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ABSTRACT 

We have identified 469 Mgll AA2796,2803 doublet systems having W r > 0.02 A in 252 Keck/HIRES and 
UVES/VLT quasar spectra over the redshift range 0.1 < z < 2.6. Using the largest sample yet of 188 weak 
Mgll systems (0.02 A <W r < 0.3 A), we calculate their absorber redshift path density, dN/dz- We find clear 
evidence of evolution, with dN/dz peaking at z ~ 1 -2, and that the product of the absorber number density and 
cross section decreases linearly with increasing redshift; weak Mgll absorbers seem to vanish above z — 2.7. 
If the absorbers are ionized by the UV background, we estimate number densities of 10 6 - 10 9 per Mpc 3 for 
spherical geometries and 10 2 - 10 5 per Mpc 3 for more sheetlike geometries. We also find that dN/dz toward 
intrinsically faint versus bright quasars differs significantly for weak and strong (W r > 1.0 A) absorbers. For 
weak absorption, dN/dz toward bright quasars is ^25% higher than toward faint quasars (10 a at low redshift, 
0.4 < z < 1.4, and 4 a at high redshift, 1.4 < z < 2.34). For strong absorption the trend reverses, with dN/dz 
toward faint quasars being ~ 20% higher than toward bright quasars (also 10 a at low redshift and 4 a at 
high redshift). We explore scenarios in which beam size is proportional to quasar luminosity and varies with 
absorber and quasar redshifts. These do not explain dN/dz's dependence on quasar luminosity. 

Subject headings: quasars: absorption lines 



1. INTRODUCTION 

Quasar absorption line systems are an extremely useful 
means of statistically constraining the various scenarios of 
metal enrichment, inflow and outflow, ionization conditions, 
kinematics, and gas structure within galaxies and the IGM. 
Various Mgll absorption line studies have concluded that 
these s ystems are cosmolog i cally distributed (lLanzetta et all 
119871; [Sargent et all Il988i ISteidel & Sargent) 1 19921) . and 
numerous subsequent studies have identif ied specific galaxies 
associated with Mgll absorption (s e e iBergeron & Boissel 



1991 



1997 



ISteidel. Dickinson. & Perssonl Il994t ISteidel et al.l 
iGuillemin & Bergeronl 1 1997F IChurchill et all 120051: 
Chen & Tinkerll2008l: iKacprzak et al1l2008l: iBarton & Cookej 
2009 1 ; IChen et alJl2010l:lKacprzak et al.ll201 U IChurchill et all 
2012tlNielsenetalJl2012l) . 

For the following discussion, we adopt the terms "weak", 
"intermediate", and "strong" to refer to absorbers having .02 
A < W r < 0.3 A, 0.3 A < W r < 1.0 A, and W r > 1.0 A, re- 
spectively, where W r is the rest frame equivalent width of 
the Mg II A2796 transition. Mg Il-selected gas probes a wide 
range of Hi column density environments. Weak Mgll ab- 
sorption in particular, which samples optically thin gas over a 
large span of cosmic time, has been prop osed to sample dwarf 
or LSB galaxies a s well as the IGM Churchill etaLI 119991; 
Rigbv etaT] |2002|) . However, there are instances in which 
weak Mg II is identifi ed in the circumgalact i c medium of "nor- 
mal" b right galaxies (Churchil l et al.l20T2l) . iMilutinovi c et all 
(2006), using an ionization model, concluded that filamentary 
and sheetlike IGM structures host at least a porti on of weak 
Mg II absorption. Moreover, Niels en et al.l (1201 2l) argue that 
W r < 0.1 A absorbers likely reside in the IGM. Photoioniza- 
tion modeling of weak M g II systems with associated Civ led 
iLvnch & CharltorJ ((2007) to argue for a scenario of a shell ge- 
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ometry as might be expected for supernova remnants or high 
velocity clouds moving in a hot corona. 

Better insight into the nature of the structures selected 
by weak Mgll absorbers remains elusive and partially mo- 
tivates this study. Measurements of the redshift path den- 
sity dN/dz of weak Mgll absorbers place important con- 
straints on n(z)<r(z ), the product of their number density 
and cross section (IChurchill et al.l fl999l iRigby et al.l l2002t 
iNaravanarT et al. 2007|)1 In the same vein, probing discrep- 
ancies in redshift path densities of Mg II absorption in vari- 
ous equivalent width ranges based on differences in the back- 
ground sources may reveal information about these absorb- 
ing structures. Despite the evidence indicating that absorption 
line systems are primarily associated with intervening gas as 
opposed to the background source itself, several studies have 
nevertheless revealed major discrepancies, depending on the 
nature of the background source, in the incidence of these ab- 
s qrbers per unit red s hift. 

IStocke & Rector! (1 19971) . in a study of strong Mgll ab- 
sorbers in BL Lac objects, observed a redshift path density 
of 4-5 times greater than that expected from quasar surveys. 
iProchter eT al. (2006b) compared quasar and GRB sightlines 
and found the latter to have a factor of ~ 4 excess in the red- 
shift path density o f strong absorbers. A similar study by 
Verga niet alJ d2009) found a factor of ~ 2 excess of strong ab- 
sorbers, but no excess of interm ediate absorbers, in the GRB 
sample. Bergeron et al. (2011) reported a factor of ~ 2 ex- 
cess of both strong and intermediate absorbers in blazar ver- 
sus quasar sightlines. It should be noted, however, that in 
that study the authors did find marginal statistical evidence 
in the strong sample of an excess occurring nearer the blazar 
(even after excluding absorbers having velocity separations 
of less than 5,000 km s" 1 from the emission redshift). Finally, 
in a study analyzin g Mgll ab s orbers down to W r = 0.07 A 
in GRB sightlines, Tejos et al J (120091) found a factor of ~ 3 
overabundance of strong absorbers and a ~ 30% reduction 
for 0.07 < W r < 1.0 A compared to studies of quasar sight- 
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lines. The latter discrepancy, however, was deemed insignif- 
icant s ince the results w ere consistent at the 1 a confidence 
level. iTeios et all d2009l) have so far presented the only study 
to compare redshift path density differences in weak Mg II ab- 
sorbers based on different types of background sources. 

In this paper we present the largest study of weak Mgll 
absorption to date, and for the first time provide a parameter- 
ized fit to the redshift path density evolution. We also present 
our findings of differential absorber redshift path density, 
based on the absolute magnitude of the background quasar, 
for weak, intermediate, and strong Mgll systems. In §|2]we 
present our data, and in §[3]our results. We discuss plausible 
interpretations in § @] and conclude in § [5] The cosmological 
parameters Hq = 70 km s" 1 Mpc" 1 , Q, m = 0.3, and J7a = 0.7 are 
adopted throughout. 

2. DATA AND SUBSAMPLES 

We have searched 252 Keck/HIRES and UVES/VLT 
quasar spectra for Mgll AA2796,2803 doublet absorp- 
tion. All s ystems were objective ly i dentified using 
the methods of [Schneider etail (11993b and IChurchill et al.1 
d!999h. Further details are provided in Evans| feOllF and 
lEvans. Churchill & Murphvl (|2012|) . The search space omit- 
ted redshifts within 5,000 km s -1 of the quasar emission red- 
shift, Zem, and blueward of the Ly a emission of the quasar. A 
total of 422 absorbers comprise our sample. 

We divided the absorber s into three subsampl es using his- 
tori cally motivated weak (|Churchill et al.l [1999), intermedi- 
ate dSteidel & Sargent! 1 19921; iNestor et al.ll2005l) and strong 
dSteidel & Sargenll992l: lNestor et al.ll2005l) equivalent width 
ranges. In order to determine whether our absorber sub- 
samples are consistent with being cosmologically distributed 
along the lines of sight to the quasars, we performed the 
second test of Bah call & Peebles! d!969t) . Performing the 
Kolmogorov-Smirnov (KS) test for each of our equivalent 
width ranges, we could not rule out that their distributions are 
consistent with being cos mological 4 . This result is in agree- 
ment with earlier studie s (Lanzetta et aLl ll987l: [Sargent et al. 
[19881: ISteidel & Sargentlll992l) . 

Since none of the quasars in our sample were observed with 
a priori knowledge of weak absorption, our survey is unbi- 
ased for this population. However, since the quasar sample 
is drawn from a broad range of targeted science programs, 
there is the possibility of bias in the intermediate and strong 
absorb ers, which are known to sometimes be associated with 
DLAs dRao & T urnshek 2000), or which in some cases were 
already known due to previous lower resolution surveys. 

To examine whether our intermediate and strong absorber 
samples are consistent with an unbiased population, we com- 
pared our measured equivalent width distributions, f(W r , 
w here W* is the chara cteristic W r , to the distribution measured 
bv lNestor et al.l d2005l). Using the KS t est for the three redshift 
bins measured by iNestor" et al. (2005), we obtained P(KS) = 
0.053 (0.36 < z < 0.87), 0.621 (0.87 < z < 1.31) and 0.322 
(1.31 < z < 2.27), respectively. For the full redshift range en- 

3 http://astronomy.nmsu.edu/jlevans/phd 

4 We calculated the Bahcall & Peebles 1 1969) Y parameter of each Mgll 
system, where < Y < 1, with 7 = representing a doublet at the min- 
imum observed redshift included in the search, and 7 = 1 representing a 
doublet at the ma ximum observed redshift. Following the formalism of 
Steidel & Sargent 1 1992), the sensitivity function g(Y) of the survey was then 
calculated. This is a measure of the number of lines of sight in the survey 
in which a system of a given minimum W r could have been detected at each 
value of Y. The Y distributions of each absorber sample were then statistically 
compared to g(Y). 



TABLE 1 

dN/dz for Mgn Surveys of 0.02 A <W, < 0.3 A 



Survey 


0.4<z<0.7 


0.7<z<1.0 


1.0<z< 1.4 


1.4<z<2.4 


CRCV99" 


1.43±0.21 


1.84 ±0.26 


2.19±0.80 




NMCK07 b 


1.06 ±0.10 


1.51 ±0.09 


1.76 ±0.08 


1.06 ±0.04 


this survey 


0.74 ±0.02 


1.08 ±0.02 


0.95 ±0.02 


0.67 ±0.01 



a IChurchilletaI]fT999l) 
b lNaravanan et alji2007l) 



compassing all three bins, we obtained P(KS) = 0.1 15. Even 
in the case of the lowest value of P(KS), corresponding to the 
lowest redshift range, the two populations are not inconsistent 
with each other to even a 2 a level. We thus proceed under the 
assumption that our sample of absorbers is a fair sample. 

For our analysis in §|3] we obtained absolute B-band and ap- 
parent magnitudes (primaril y B, V, and R) of the quasars. The 
majority were obtained from lVeron-Cetty & V eron (200J}; 16 
were obtained from the NASA/IPAC Extragalactic Database; 
and for two quasars, the magnitudes could not be determined 
so these lines of sight and their absorbers were omitted from 
analysis for which these quantities were required. 

3. RESULTS 

Following the formalism of lLanzetta et al.l dl 987). modified 
to account for the doublet ratio (Churchill et al . 1999), we cal- 
culated the number of absorbers per unit redshift, dN/dz. 

3.1. Weak Absorber Redshift Path Density 

For our full redshift range, 0. 1 < z a bs < 2.6, the cumulative 
redshift path is AZ = 231 and is ~ 100% complete to a 5 a 
equivalent width sensitivity of W,-(2796) = 0.05 A and ~ 82% 
complete to a 5 a equivalent width sensitivity of W, (2796) = 
0.02 A. Over the redshift range 0.4 < z a bs < 2.4, the extent 
of a study by Nara yanan et al.l (120071) . we have AZ = 213, 
compared to their 70; and for 0.4 < z a bs < 1 -4, the extent of 
a study by IChurchill et al.l d!999l) . we have AZ = 148, com- 
pared to their 17. These represent the largest two previous 
weak Mgll surveys. Our larger cumulative redshift path re- 
flects the larger number of lines of sight included in our sur- 
vey. IChurchill et al.l (1 19991) surveyed 26 HIRES quasar spec- 
tra and found 30 weak systems, while Naraya nan et al.l (120071) 
surveyed 81 UVES quasar spectra and found 112 weak sys- 
tems. In our survey we identified 188 weak systems. 

We calculated dN/dz for weak systems in four redshift 
bins in order to facili tate comparison with the works of 
iNarayanan eTall d2007l) and IChurchill et al.l (1 19991) : the bins 
and results are shown in Table [TJ 

All three studies obtained different results, with the trend 
being that the dN/dz values have decreased with larger sur- 
vey size. All three used the identical code (SEARCH, 
IChurchill et al.l [l999h for line and candidate identification, 
and all three ostensibly used the same algorithms in deter- 
mining the redshift paths for each system. To test for possi- 
ble differences in the calculations 5 , we ran our code used for 
thi s survey on the identi cal spectra and set of systems used 
by IChurchill et all (1 19991) and compared their AZ and dN/dz 
results with those of the r eproduced study. The result was 
that the AZ values from the lChurchill et al.l dT999) study were 
~ 30% larger than our reproduced values from their data, and 

5 The Churchill et al. (1999) spectra were a subset of the quasars we 
searched. 
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the dN/dz values from the original study were thereby low- 
ered compared to our reproduced study. This suggests that our 
redshift path ca l culat ions are more conservative than those of 
IChurchill et aD (fl999l If our redshift paths had been calcu- 
lated exactly as theirs, the dN/dz result of this survey would 
presumably have been lower, further widening the discrep- 
an cy with previous works. A similar duplication of the survey 
of Naraya nan et alj (120071) could not be performed because 
a significant number of the authors' quasar spectra were un- 
available to us. 

In an attempt to find differences in the quasar samples that 
could possibly lead to the discrepant dN/dz results among 
the three works, we investigated the quasar apparent magni- 
tude distributions. The KS test was performed between all 
possible pairs of the th ree surveys. Our survey and that of 
iNarava nan et alj (120071) exhibit remarkably similar distribu- 
tions; it could not be ruled out to greater than a 1 a confidence 
level that their apparent magnitudes had been drawn from the 
same popul ation, and their median values were both 17.5. In 
contrast, the Church ill et al.l ([1999) survey differed from each 
of the other two to a confidence level of 4a and had a median 
apparent magnitude of 16.3. A uniform set of apparent mag- 
nitudes in the same band was not available, however, mak- 
ing the se comparisons uncertain. The smaller IChurchill et al.l 
(1999) survey was undoubtedly overall biased toward brighter 
quasars. 

We likewise investigated the quasar absolute Z?-band mag- 
nitude distributions, but when each survey was tested against 
the other two, it could not be ruled out to even a 1 a confidence 
level that their distributions had been drawn from the same 
population. Similarly, no substantial differences were found 
in the overall median quasar absolute magnitude of the thre e 
surveys: -28.5,-28.7 and -28.6 for IChurchill etal.l (119991) . 
INaravanan et alj d2007). and this study, respectively. 

The reasons behind the different dN/dz results among the 
three works remain unclear; future studies will hopefully re- 
solve the weak Mgll puzzle. It may be that the inclusion or 
exclusion of weak systems very close to the limiting equiva- 
lent width may play a role (Anand Narayanan, private com- 
munication), as well as differences in the codes used, since 
the calculation of AZ is extremely sensitive for the weakest 
systems. 

3.2. Redshift Evolution 

The number of absorbers per unit redshift is the product 
of the proper number density of absorbers and their proper 
geometric cross section. In the standard cosmological model, 
the no-evolution expectation (NEE) for the redshift number 
density can be parameterized as: 



dN 
~dz~ 



c (1+z) 2 

Ho ^/n m (i+z) 3 +n A 



(i) 



where «o is the mean comoving number density of absorbers 
and do is the mean comoving geometric absorber cross sec- 
tion. As shown in Figure QJ, our observed dN/dz departs 
strikingly from the NEE (dashed curve), which was normal- 
ized to the mean Zabs, 1-098, and to the overall dN/dz, 0.83, 
of our weak sample over the redshift range 0.1 < tabs < 
2.6. The general beha vior of dN/dz is in agreement with 
INaravanan et alj (120071) in that it peaks between 1 .0 < z < 14 
and then decreases toward higher redshift. 

If the product na varies as a function of redshift, dN/dz 
may depart from the no-evolution expectation. The quantity 



n(z)a(z) can be written as 

n(z)a(z) = n a f(z), 



(2) 



where f(z) is a nonnegative function that parameterizes the 
evolution of dN/dz. 

Figure Q}? depicts our weak dN/dz result divided by the 
NEE. The data clearly motivate a linear fit; this was achieved 
using a function of the form 



f(z)=l-a(z-z*), 



(3) 



where a is the slope and z* is the function normalization. The 
result, a = 0.69 ±0.02 and z* = 1.29 ±0.05, is shown as a 
solid line in Figure 03?. 




FIG. 1. — (a) shows the weak dN/dz (open circles) as well as the NEE 
(dashed line). — (b) shows the same data divided by the NEE. The solid line 
is the best linear fit to this ratio. In (a), this fit was then multiplied by the 
NEE and plotted (solid curve). 

Finally, the fit to f(z) was multiplied by the NEE, yielding 
a fit to the dN/dz result in FigureQJ (solid curve) that can be 
expressed as 

(l+z) 2 [l-a(z-z*)] 



dN c 

~T = TT n ° ao ~ ; 



(z<2.74). (4) 



Weak Mgll absorbers seem to vanish at high redshift; this 
result is discussed in §|4] 

3.3. Differential Absorber Redshift Path Densities by 
Absolute Magnitude 

Motivated by studies that found differences in dN/dz re- 
sults based on background object as discussed in § Q] as well 
as our finding of differing apparent magnitude distributions 
among the weak Mg II studies as discussed in § 13.11 we at- 
tempted to discern some intrinsic difference that might affect 
the observed dN/dz- Using the absolute B-band quasar mag- 
nitudes of our survey, which had a range of -32.2 < Mb < 
-19.8, we divided our quasars into "bright" and "faint" sub- 
samples according to the median, (Mb) = -28.6. The bright 
subsample has a median of -29 .4, and the faint subsample has 
a median of -27.4. 

Figure |2] shows our weak dN/dz results for the bright and 
faint samples binned as in Figure Q] The redshift path density 
of the bright sample is clearly higher than that of the faint 
sample at all redshifts. The lowest redshift bin of the weak 
bright sample contained only five systems and we believe this 
caused the dN/dz result in that bin to be less reliable. 

Figures |3]2-yt plot our dN/dz results for weak, interme- 
diate and strong Mgll absorption. The bins (0.4 < Zabs < 
1.4 and 1.4 < z a bs < 2.34, "low" an d "high" redshift) were 
selected based on previous studies ([N arayanan et al.l 120071: 
iNestor et al.l 12001 IChurchill et al.1 [19991 Similarly, Fig- 
ures [3j^43f show dN/dz calculated for the faint and bright 
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TABLE 2 
dN/dz by Quasar Luminosity 



simple 


weak 


sig 


intermediate 










0.02 < W r < 0.3 


lev 


0.3 < W r < 1.0 




W,- > 1.0 


lev 




0.4 < z abs < 1.4 


all 


0.931 ±0.006 




0.715 ±0.005 




0.449 ± 0.003 




hi i trlu 


1.040 ±0.016 




0.692 ± 0.010 




0.406 ± 0.006 




faint 


0.836 ±0.011 




0.730 ± 0.009 




0.481 ±0.006 




bright/all 


1 .1 17 ± 0.019 


6.2<t 


0.967 ±0.016 


2.1 (7 


0.905 ±0.015 


6.3 a 


faint/all 


0.897 ±0.013 


7.9 er 


1.021 ±0.014 


1.5(7 


1.072 ±0.015 




faint/bright 


0.803 ± 0.020 


9.9(7 


1.056 ±0.019 


2.9(7 


1.185 ±0.019 


9.7 cr 


I* < Zabs <234 


all 


0.686 ±0.011 




0.696 ±0.011 




0.631 ±0.010 




brighl 


0.732 ±0.018 




0.730 ±0.018 




0.584 ±0.014 




faint 


0.579 ± 0.028 




0.619 ± 0.030 




0.714 ±0.034 




brighi/all 


1.068 ±0.031 


2.2ct 


l.O5O±0.O31 


1.6(7 


0.926 ± 0.027 


2.7 a- 


faint/all 


0.845 ± 0.043 


3.6(7 


0.890 ± 0.045 


2.4(7 


1.131 ±0.057 


2.3 cr 


faint/bright 


0.791 ±0.054 


3.9(7 


0.848 ± 0.054 


2.8(7 


1.222 ±0.053 


4.1 (7 



intermediate 





FIG. 2. — dN/dz of the weak sample for the bright (red triangles), faint 
(blue squares), and combined (dotted crosses) samples. In the bright sample, 
the lowest redshift data point has a value of 1 .53 ± 0.47. 

quasar samples. Figures |3^-[3}' show the ratio of the dN/dz 
results of the faint quasar sample to that of the bright. These 
values are listed in Table [2] for each W r and redshif t range . 
Though the weak dN/ dz results of IChurchill etafl (119991) . 
iNarayanan et al.1 (120071) . and this study all differed as dis- 
cussed in § 13.11 the relative values among our own faint, 
bright and all quasar subsamples are robust, having been cal- 
culated in a consistent manner. Similarly, though our interme- 
diate and strong sam ples have W r distributions consistent with 
iNestor et al.l (120051) as mentioned in § [21 our dN/dz results 
for these samples are high er than those of Nestor et al. (2005) 
and lLundgren et all (l2009h . This is expected since some of the 
quasars in our survey were targeted for their known W r > 0.3 
A Mgll absorption; however, the dN/dz ratios between our 
magnitude bins are robust. 
For weak systems, the dN/dz values toward bright quasars 



FIG. 3. — (a-c) show dN/dz of the weak, intermediate and strong samples, 
respectively, in two redshift bins (0.4 < z a bs < 1.4 and 1.4 < z a bs < 2.34). 
— (d-f) show dN/dz for the bright (red triangles) and faint (blue squares) 
subsamples. The points were left open to better display the error bars. — 
(g-i) show the faint to bright dN/dz ratio. 

are ~ 25% higher than for the faint in both redshift ranges, 
while for strong systems, dN/dz is ~ 20% higher toward faint 
quasars than toward bright in both redshift ranges. The faint to 
bright dN/dz ratio, (dN/dz)//(dN/dz)b, departs for the weak 
systems from unity at the ~ 10 a level for low redshift and at 
the ~ 4 a level for high redshift. The strong absorber ratios 
are similarly significant (see Table |2}, though in that case the 
faint quasar dN/dz values are higher, rather than lower, than 
those of the bright sample. In the case of the intermediate 
absorbers the ratio is consistent with unity within 3 a. 

In the calculation of dN/dz the sensitivity of each sp ectrum 
is accounted for (see Lanzetta et al. 1987; Evans 201 1), elimi- 
nating the possibility that the higher values of the weak bright 
sample compared to the weak faint sample might result from 
higher signal-to-noise ratios. 

4. DISCUSSION 

Based on our weak dN/dz result (see Figure [T), the cos- 
mic number density, geometric cross section, or both, of weak 
Mgll absorbers appear to be evolving. The apparent dropoff 
in our fit toward z = may be overly steep; it is possible that 
the weak dN/dz peaks at z ~ 1, declines slightly and then 
levels off toward the present. However, it is a first attempt to 
characterize weak Mg II absorber evolution using a functional 
form. Our result predicts that no such absorbers exist above 
z ^2.7. 

INarayanan et al.l (120071) speculated that the apparent 
paucity of weak Mg II above z ~ 2 might be due to the high 
redshift analogs of low redshift weak Mg II absorption being 
associated with strong Mg II. In this scenario, weak Mg II ab- 
sorption at high redshift would be in the kinematic vicinity 
of strong Mg II and thus would not be recognized as isolated 
weak absorption. However, we have compared the high ve- 
locity weak kinematic subsystems o f strong Mg II subsystems 
to isolated weak Mgll absorption dEvansll201 ll) . Morpho- 
logically these two types of profiles often appear very simi- 
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lar, but a KS test of their rest equivalent width distributions 
revealed that they are actually two distinct populations to a 
99.98% confidence level, or greater than 3 a. A KS test of the 
distributions of flux decrement-weig hted velocity spreads, lu v 
(IChurchill& Vogtll200lL lEvansll201 ll) indicated to a greater 
than 6 a confidence level that the two populations are unique. 

The dN/dz evolution we detect in weak Mgll may be due 
to changes in gas structure or ionization conditions; neither 
we nor other studies find the same falloff in intermediate 
and strong Mg H absorption (|Nestor et al.l2005t iProchter et al.l 
l2006at iLundgren et al.l 12009b up to our maximum redshift 
of 2.6. Matejek & Simcod d2012h do report a decline in the 
strong population above z ~ 3, and note that this peak corre- 
sponds to that of the SFR. Our weak dN/dz result, which ex- 
hibits no such peak, may provide indirect evidence that a sub- 
stantial fraction of these absorbers resides in the IGM, since 
their evolution appears not to correspond to star formation. 

Using Cl o udy 08.00 photoionization modeling 
dFerland et al.1 [1998), we investigated the evolution of 
weak Mg II absorber sizes R(z ) and cosmic number densities 
n(z) (for additional details see lEvansll2011l) . In this scenario 
Mgll selects relatively dense cloudlets embedded within 
plane parallel slabs of gas. We modeled optically thin clouds 
having a range of hydrogen number densiti es based on past 
Mgll photoionizatio n modeling results (iRigbv et al.l 120021; 
Berge ron et al.l l2002). We assumed an ultraviolet background 
model that varies as a f unction of z, following the work of 
lHaardt & Madaul d!996l) . which includes the contribution of 
galaxies. Though we examined a grid of clouds with a range 
of Hi column densities and metallicities, we discuss here 
clouds having A/(Hl) of 10 16 cm" 2 and a metallicity of 0.1 
solar. For weak Mgll, logjV(Hl) is constrained to the r ange 
15.5-17.0 cm" 2 (IChurchill et alJ2000t IPJgbvltini2002h . 

The resulting cloud thicknesses, which we interpreted as 
absorber sizes R(z) and which are governed by the ionizing 
background, peak at z ~ 2 and then decline toward the present, 
as shown in Figure HJj. Assuming spherical clouds, the cor- 
responding absorber cross sections a(z), combined with our 
weak dN/dz constraint using the fit of Equation [4] trans- 
late into cosmic absorber number densities n(z) that increase 
monotonically toward the present (Figure HI?). The absorber 
sizes produced by this model are on the order of a parsec, 
and yield absorber number densities on the order of 10 6 - 10 9 
Mpc" 3 , for the middle range of values. T his corresponds to 
10^-10 12 absorbers per L * galaxy for z < IdFaber et al.l2007l) 
as we ll as for 1 < z < 3 (Red dv & Steidell2009l; lOesch et alj 
l20Toh . 

If the clouds are not spherical, but instead the transverse 
extent Rj scales with cloud thickness according to a factor 
/3 such that Rj = /3R(z), then n{z) would scale as /3~ 2 . For 
(3 = 100, n(z) would then be reduced by a factor of 10 4 , which 
yields 10 2 - 10 s weak absorbers per Mpc 3 . It should be noted 
that changing the model's A/(Hl) would change R(z) in direct 
proportion, while n(z), using our dN/dz constraints, would 
vary as N(Hi)~ 2 . 

Our Cloudy model is suggestive of a condensation mech- 
anism into sheet or filament structures characteristic of the 
IGM. Sheetlike geometries require far fewer weak absorbers 
than do spherical geometries per L* galaxy, and therefore we 
consider these to be a more realistic scenario. Although we 
do not fully explain the nature of weak Mg II absorption, this 
exercise does provide a limiting case in which to couch the 
phenomenon. Our absorber size estimate is not far from that 
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FIG. 4. — (a) plots the evolution of Cloudy absorber size, log i?(z). (b) 
plots the number density of absorbers, log n(z). The dashed curves represent 
log uh = —0.5 cm -3 ; the dotted curves, log = —1.0 cm~ 3 ; the solid curves, 
log riH = —1.5 cm -3 ; the long dashed curves, log iih = -2.0 cnT 3 ; and the 
dashed-dotted curves, log nu = — 2.5 cm -3 . 

of Rig bv et ail d2002l) . who concluded that a sheetlike struc- 
ture containing many embedded ~ 10 pc absorbers was re- 
quired to account for t he o bserved dN / dz . Rece nt findings by 
IChurchill etall d2012l) and lNielsen et all d2012l) indicate that 
0.1 < W r < 0.3 A, but not 0.02 < W r < 0.1 A, absorbers are 
found in the circumgalactic medium of normal galaxies at im- 
pact parameters of less than 200 kpc. They suggest that the 
weaker of these two populations may then reside primarily in 
the IGM. 

Our dN/dz results toward faint versus bright absolute mag- 
nitude quasars reveal that the weak absorbers have a higher 
redshift path density in the bright quasar sample than in the 
faint at both low and high redshifts. For the strong absorbers, 
the opposite is true. The intermediate absorbers appear to fol- 
low no clear trend, and may represent an equivalent width 
range where the effects leading to the weak and strong differ- 
entials mostly cancel. 

Following the discovery by IProchter et al.l (l2006bl) of the 
quasar-GRB strong Mg II discrepancy, various researchers 
have attempted to explain the phenomenon. Fran k et all 
d2007l) modeled the effects of Mg II absorber size and impact 
parameter on observed equivalent width and concluded that 
the dN/dz discrepancy may be due to the larger beam sizes 
of quasars versus GRBs, the latter of which they state are on 
the order of the sizes of cloud cores. The authors also pre- 
dicted that different luminosity populations of quasars should 
contain different incidences and strengths of intervening ab- 
sorbers; 

iPorciani et al.1 d2007l) countered this differential beam size 
argument, noting that no unsaturated Mg II doublets have been 
observed having a doublet ratio of one, as would be expected 
in the case of partial covering of a quasar beam by an ab- 
sorber. They stated that magnification bias could explain the 
discrepancy, and that dust obscuration bias and association of 
absorbers with the circumburst environment could also par- 
tially account for it. 

A s tatistical study was also conducted by iPontzen et alj 
d2007l) to look for systematically lowered Mgll equivalent 
widths over quasar broad line emission regions, which are 
substantially larger than quasar continuum regions; no sig- 
ni ficant difference was fo und. 

ICucchiara et al.l d2009l) cite an intrinsic origin as a possible 
explanation for the GRB excess. Mgll absorbing gas could 
be ejected at relativistic velocities and masquerade as an in- 
tervening absorber; however, the authors note that the pres- 
ence of Mg I absorption in these systems, as well as the lack 
of fine structure transitions that are expected in the vicinity of 
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a GRB, cast doubt on this theory. Verga niet al.l d2009) concur 
that the excess could be intrinsic, and estimate required ejec- 
tion velocities of 10, 000-25, 000 km s . They also consider 
gravitational lensing to be a viable mechanism to account for 
the discrepancy. 

In two stud ies of Civ absorbers towar d GRBs, 

ISudilovskv et al.l d2007l) and iTeios et al.l d2007l) reported 
no exc ess incidence over quasar sightlines. ISudilovs kv et al.l 
(120071) speculated that the difference in the cases of Mgll 
versus Civ absorbers arose partially because the former 
introduced more dust extinction than the latter. However, 
dust extinction in Mg ll absorbers was subsequently modeled 
(ISudilovskv et al.ll2009t) . and the authors concluded that the 
effect could only account for ~ 10% of the quasar-GRB 
di screpancy. 

Tejo s et aTl (|2009) rejected an intrinsic origin for excess 
GRB Mg II absorbers due to the lack of both excess C IV ab- 
sorption and excess weak and intermediate Mg II absorption, 
and ins tead favored gravita tional lensing as the relevant mech- 
anism. IWyithe et al.l (1201 ll) modeled gravitational lensing in 
quasar and GRB sightlines and concluded that it was a fea- 
sible explanation for the excess, but that further GRB data 
were necessary to support or refute their findings. They noted 
that afterglows in which strong Mgll systems are found are 
brighter than average, implying a greater lensing rate. 

Through th eir modeling of extinc ti on cu rves toward quasars 
and GRBs, iBudzynski & Hewettl ( 1201 lh determined that 
dN/dz toward quasars would be significantly higher if cor- 
rected for dust, and that the correction varies with redshift. 
The discrepancy compared to GRBs arises, the authors state, 
because their absorber redshift distribution is shifted higher 
than that of quasars, resulting in less loss of detected ab- 
sorbers. They calculated that this effect could account for a 
factor of two excess in the GRB dN/dz. 

Keeping these previous studies of the quasar-GRB discrep- 
ancy in mind, and in an attempt to understand the possibly 
related phenomenon we have uncovered, we attempted to find 
some other metric within our faint and bright quasar samples 
that would shed light on these issues. The redshift path den- 
sity depends on the integrated equivalent width distribution, 



dN 
dz 



f(W,W*)dW. 



(5) 



We therefore compared equivalent width distributions for our 
faint and bright quasar samples for the various W r ranges as 
well as for low, high, and all redshifts. 

We also studied equivalent width distributions binned by 
relative beam sizes using Mb as a proxy, i.e. assuming that 
the square of the source radius R s is proportional to the B- 
band luminosity of the quasar. The ratio of the source radius 
for quasar i relative to the median radius for the full quasar 
sample can then be written (Shak ura & Sunv aev 1973;) as 



R 



"V _ 1 n-0.2(M,-(M if )) 



(6) 



We then studied the effect of changing beam size with redshift 
due to cosmology. We calculated the ratio of the relative beam 
size of quasar i at the redshift of absorber j to the cross section 
of the source: 



Vb(Zj) 



D A (z s ,i) 



(7) 



where a s j = nR 2 si , Da{Zj) is the angular diameter distance at 
the absorption redshift of system j, and Da(z s j) is the angular 



diameter distance at the source redshift of quasar i. Finally, 
we examined the combined effect of source size and cosmol- 
ogy, by using Equation|6]to scale er^. 

Using the KS test, none of these equivalent width distribu- 
tions yielded significant differences (of at least 3 a) between 
the faint and bright quasar samples. Since the dN/dz discrep- 
ancy is in this case of a smaller magnitude than in the case of 
the quasar-GRB phenomenon, it may require a larger data set 
to discern the reaso ns behind the observations . 

The findings of Budzynski & Hewettl d201 lh do offer an in- 
triguing possibility by relating absorption redshift distribu- 
tions to dN/dz. Our faint absorber samples do have lower 
median Zabs values than our bright samples across all three 
equivalent width ranges, probably a result of the correlation 
whereby intrinsically more luminous quasars tend to be se- 
lected at higher redshifts. This result only supports the au- 
thors' dust argument in the case of our weak absorbers, the 
only sample in which the absorption incidence is significantly 
higher toward bright quasars than faint. For our weak sample, 
the median absorption redshift is 0.87 in the faint subsample 
and 1.17 in the bright. KS testing, however, revealed that it 
could not be ruled out to a greater than 98.32% confidence 
level that the faint and bright subsamples are drawn from the 
same underlying z a bs distribution. 

Though several authors have argued for an intrinsic ori- 
gin for excess strong Mgll absorption toward GRBs versus 
quasars, this does not appear to explain the dN/dz discrep- 
ancy in_flie_£ase_of_our bright and faint quasar populations. 
Our lBahcall & Peebles! (|T969) testing revealed absorber dis- 
tributions consistent with cosmological within all equivalent 
width and quasar absolute magnitude subsamples, as well as 
in the aggregate populations. The velocities of Mg Il-selected 
gas ejected from a quasar would have to reach large fractions 
of the speed of light in order to pass for intervening systems. 
It therefore seems highly unlikely that significant intrinsic ab- 
sorption could be present in our sample. 

5. CONCLUSION 

We have found in a survey of 252 quasar spectra that the 
incidence of weak Mg II absorption evolves markedly, that it 
peaks at z ~ 1 -2, and that it is fit by a function that is a product 
of the no-evolution expectation with a linear function. Our 
linear fit to the ratio of our dN/dz data to the NEE resulted in 
a slope of a = 0.69 ± 0.02 and a normalization of z* = 1 .29 ± 
0.05 for the function f(z) = l-a(z-z*). Our dN/dz result 
predicts that no weak Mgll absorbers exist above z — 2.7. 

We find that when our quasar survey is segregated by abso- 
lute magnitude, weak Mgll dN/dz is significantly lower in 
the faint subsample than in the bright, with faint to bright 
dN/dz ratios of 0.80 ± 0.02 at low redshift and 0.79 ± 0.05 
at high redshift. In contrast, strong Mgll dN/dz is signifi- 
cantly higher in the faint subsample than in the bright, with 
faint to bright dN/ dz ratios of 1 . 19 ± 0.02 at low redshift and 
1 .22 ±0.05 at high redshift. Intermediate equivalent width ab- 
sorbers exhibited dN/dz ratios consistent with unity within 3 
a. At this time it is uncertain whether these results stem from 
some intrinsic property of the quasars, from some difference 
in the intervening sightlines, or from some combination of 
these factors. 
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